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Abstract 

We consider the fully constrained version of the next-to-minimal supersymmetric 
extension of the standard model (cNMSSM) in which a singlet Higgs superfield 
is added to the two doublets that are present in the minimal extension (MSSM). 
Assuming universal boundary conditions at a high scale for the soft supersymmetry- 
breaking gaugino, sfermion and Higgs mass parameters as well as for the trilinear 
interactions, we find that the model is more constrained than the celebrated minimal 
supergravity model. The phenomenologically viable region in the parameter space of 
the cNMSSM corresponds to a small value for the universal scalar mass mo: in this 
case, one single input parameter is sufficient to describe the phenomenology of the 
model once the available constraints from collider data and cosmology are imposed. 
We present the particle spectrum of this very predictive model and discuss how it 
can be distinguished from the MSSM. 
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The next-to-minimal supersymmetric standard model (NMSSM) [1-3], in which the 
spectrum of the minimal extension (MSSM) is extended by one singlet superfield, was 
among the first supersymmetric (SUSY) models based on supergravity-induced SUSY- 
breaking terms. It has gained a renewed interest in the last decade, since it solves in 
a natural and elegant way the so-called fi problem [4] of the MSSM; in the NMSSM 
this parameter is linked to the vacuum expectation value (vev) of the singlet Higgs field, 
generating a value close to the SUSY-breaking scale. Furthermore, it leads to an interesting 
phenomenology as the MSSM spectrum is extended to include an additional CP-even and 
CP-odd Higgs state as well as a fifth neutralino, the singlino. 

In contrast to the non- or partially constrained versions of the NMSSM that have been 
intensively studied in the recent years [5], and which involve many free parameters, the 
constrained model (cNMSSM) has soft SUSY-breaking parameters that are universal at a 
high scale. This is motivated by schemes for SUSY-breaking that are mediated by (flavor 
blind) gravitational interactions, and leads to a more predictive model as the number of 
unknown parameters is reduced to a handful. 

In the constrained MSSM or minimal supergravity (mSUGRA) scenario [6], the uni- 
versality of the soft parameters at the grand unification scale Mgut leads to only four 
continuous parameters: the common gaugino M\i<i and scalar m masses, the trilinear 
coupling Aq, and the ratio of vevs of the two Higgs fields tan/3. (The conditions for a 
correct electroweak symmetry breaking allow to trade the two basic Lagrangian parame- 
ters /i 2 and its corresponding soft breaking term B\x by the Z boson mass Mz and tan/?, 
leaving the sign of \i undetermined.) The numerous soft parameters at low energies are 
then obtained through renormalization group (RG) evolution. The cNMSSM with univer- 
sal soft terms at the GUT scale has the same number of unknown parameters: instead 
of the cMSSM \i and B parameters, one has the trilinear couplings A and k in the Higgs 
sector. 

General features of the cNMSSM parameter space as well as its phenomenology have 
been discussed in Refs. [2,3]. These studies already revealed that the allowed range for 
the parameters Mi/2,mo and Aq is different from that of the cMSSM. Indeed while small 
values for mo are disfavored in the cMSSM as they lead to charged sleptons that are lighter 
than the neutralino x?? which is the preferred lightest SUSY particle (LSP), small tuq is 
needed in the cNMSSM to generate a non- vanishing vev of the singlet Higgs field [2]; the 
slepton LSP problem can be evaded owing to the presence of the additional singlino-like 
neutralino which, in large regions of the cNMSSM parameter space, is the true LSP. Note 
that vanishing values for m are naturally obtained in supergravity models with Kahler 
potentials of the no-scale type [7]; however, the additional no-scale prediction A = is 
difficult to realize in the cNMSSM. 

Since the early studies, bounds on the Higgs and SUSY particle spectrum from high- 
energy collider data and low-energy measurements have become more severe, while im- 
portant inputs such as the top quark mass are more accurately measured [8] . In addition, 
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tools for a more precise determination of the mass spectrum and couplings [9], and the 
cosmological dark matter relic density [10] have become available. In this letter, we re- 
investigate the parameter space of the cNMSSM in the light of these recent constraints, 
using the updated tools. 

We find a viable region of the cNMSSM parameter space very close to the no-scale 
situation: vanishing or very small mo and \A \ well below Mi/ 2 . In fact, in order to 
obtain the correct cosmological relic density for the LSP, which turns out to be a singlino- 
like neutralino [11], Aq is essentially fixed in terms of Mi/ 2 ', tan/? is determined by the 
scalar mass universality (including the additional singlet), while LEP constraints on the 
Higgs sector imply A < 10 -2 and this parameter has practically no effect on the remaining 
particle spectrum. 

Thus, remarkably, only one single parameter is sufficient to describe the phenomeno- 
logical features of the cNMSSM. This makes the model much more predictive than the 
celebrated mSUGRA model. In addition, as will be shown, the phenomenology differs 
considerably in the two scenarios. 

In this analysis, we consider the NMSSM with a scale invariant superpotential given 

by 

W = XSH u H d + ^S 3 + ... (1) 

where the two terms shown substitute the fiH u Hd term in the MSSM superpotential and 
we have omitted the usual generalization of the Yukawa interactions. The soft SUSY- 
breaking terms consist of mass parameters for the gauginos M 12 ,3, sfermions rnp LR and 
Higgs fields rriH ud and trilinear interactions Af as in the MSSM, supplemented by an 
additional scalar mass m>s and two trilinear couplings A K and A\ for the singlet field. 
Once the unification of these soft SUSY-breaking masses and trilinear couplings at the 
scale M G ut is imposed, the Lagrangian of the cNMSSM depends on the five parameters, 
Mi/2) m o, A , A and k. The correct value for Mz reduces the dimension of the parameter 
space from five to four; e.g. k can be determined in terms of the other parameters. Hence, 
the number of continuous free parameters in the cNMSSM is the same as in the cMSSM. 

For practical purposes it is convenient to adopt the following procedure: In addition 
to Mz, one allows for the five cNMSSM input parameters 

^i/2 ) m o > A) > X and tan (3 . (2) 

The parameters k and the soft singlet mass squared m| are both determined at the scale 
Msusy through the minimization equations of the scalar potential. (The vev (S) or fi e g = 
X (S) is also fixed through the third independent minimization equation, which leaves 
the sign of // e fr undetermined.) This is the procedure employed by the routine NMSPEC 
within NMSSMTools [9], which calculates the spectra of the Higgs and SUSY particles in 
the NMSSM in terms of the soft SUSY-breaking terms at Mqut (except for the parameter 
m|), tan/? at the weak scale and A at the SUSY-breaking scale Msusy- 



3 



Clearly, the soft singlet mass squared m 2 s at Mqut will in general not coincide with 
777,0 ■ However, one can confine oneself to regions in parameter space where the difference 
between m 2 s and 777.0 * s negligibly small. This condition leaves us with an effective four- 
dimensional parameter space, consistent with the considerations above. 

Let us now explore the cNMSSM parameter space defined by arbitrary values of M1/2, 
A and A, assuming fi e s positive and, for the time being, restricting ourselves to m = 0. 
For each choice of Mi/ 2 , Ao and A, tan f3 is determined by the requirement that m| at Mqut 
should be close to 777,0 = 0. In practice, we impose m|(M G uT) < (5 GeV) 2 , which typically 
requires to tune the fourth decimal of tan f3 (this should not be interpreted as a finetuning, 
since 7775 should be considered as an input parameter, whereas tan (3 is determined by the 
minimization of the effective potential). For the most relevant SM parameters we chose 
a s (M z ) = 0.1172, m b (m b ) m = 4.214 GeV and = 171.4 GeV [8]. 

For this set of parameters, we select the cNMSSM space which survives once one 
imposes theoretical requirements such as correct electroweak symmetry breaking, pertur- 
bative couplings at the high scale, the absence of tachyonic masses, a neutralino LSP, 
etc., and phenomenological constraints such as the LEP bounds on Higgs masses and 
couplings, collider bounds on the SUSY particle masses [8], experimental data from B- 
meson physics [12] and from the anomalous magnetic moment of the muon [13], and a 
relic density compatible with cosmological data [14]. For A < 10~ 2 , as it turns out to be 
the case (see below), the phenomenologically allowed region is nearly independent of the 
input A. 

Leaving aside, for the time being, the WMAP constraints on the LSP relic density 
Qh 2 , but including all the other constraints above, the phenomenologically allowed region 
in the [Mi/ 2 , A ] plane is shown in Fig. 1. 

This allowed region is bounded from below, i.e. for large absolute values of |A |, by 
the absence of a charged (generally tau slepton) LSP as in the cMSSM for mo = 0; inside 
the allowed region, the LSP is a singlino-like state. The upper bound at A ~ follows 
from the positivity of the mass squared of the singlet-like CP-odd Higgs boson, which is 
given to a good approximation by — 3kA k (S), k (S) being positive, and A K ~ A . To the 
left, i.e. towards smaller values of M1/2, the allowed region is bounded simultaneously 
by the condition that the lightest tau slepton mass must be above ~ 100 GeV from its 
non-observation at LEP, and the mass of the lightest SM-like CP-even Higgs boson above 
~ 114 GeV. To the right, i.e. towards larger values of M1/2 for which SUSY particles 
are too heavy to be produced at the LHC, no "physical" constraint appears. We believe, 
however, that this region requires more and more hidden finetuning 0. Inside this allowed 
region, tan/5 is fixed by the requirement m|(McuT) ~ 0, and turns out to be quite large 
(see Ref. [15] for earlier work on the cNMSSM at large tan f3). In Fig. 1, we have indicated 
lines corresponding to constant tan (3 = 25, 30, 33 and 35. 

lr Thc actual finetuning is hidden by our procedure which determines k in terms of other parameters 
such as Mz; this is similar to the determination of B in the cMSSM. 
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Figure 1: The viable cNMSSM region in the [Mi/ 2 , A ] plane for m ~ and A = 2 x 10 -3 , 
once theoretical, collider and cosmological constraints have been imposed. 

As a next step, we require that the WMAP constraint on the relic density of the 
Xi dark matter (DM) candidate, calculated using the program MicrOMEGAS [10], is 
satisfied. Given the actual small error bars, 0.094 < Q x ah 2 < 0.136 at the 2a level [14], 
this constraint (if satisfied at all) reduces the parameter space of any model to a lower 
dimensional hypersurface. Within the allowed region in Fig. 1, the correct relic density for 
Xi is obtained along the line close to the lower boundary (DM line) , where the mass of the 
singlino-like LSP is close to the mass of the next-to-LSP which is the lightest tau slepton 
f\, Mf ± — M x o ~ 3 to 5 GeV. Near the upper boundary, the mass of the singlino-like LSP 
becomes very small, implying a far too large relic density. For A 1 (see below), the 
singlino-like Xi LSP has very small couplings to SM particles and, thus, the main process 
which generates the correct relic density is the assisted annihilation or "co-annihilation" 
of tau sleptons, fx fx — > SM particles, which does not depend on these coupling^. 

As mentioned previously, we have checked that constraints from S-physics [12], such 
as the branching ratio for the radiative decays b — > 57, are satisfied. Moreover, we find 
that the supersymmetric contribution to the anomalous magnetic moment of the muon, 
<5a M — (g — 2) M , accounts for the ~ 3a deviation from the SM expectation [13]: along 
the DM line, 5a^ decreases from ~ 4.2 x 10~ 9 for M 1/2 ~ 400 GeV to ~ 0.2 x 10~ 9 for 
Mx/2 ~ 1-5 TeV. In view of a desired value Sa^ ~ (2.7 ±2) x 10~ 9 , the region M 1/2 < 1 TeV 
is thus preferred by this observable. 

2 Howcvcr, A values larger than ~ 10 -5 , as used here, are needed such that the hypothesis of ther- 
mal equilibrium between the LSP and the next-to-LSP near the relevant temperature can be considered 
satisfactory. 
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At this stage, one can conclude that for small mo and A, only one single parameter 
is needed to describe the phenomenology of the cNMSSM, once collider and cosmological 
constraints are imposed. Along the DM line in Fig. 1, once a value for Min is chosen, 
both the universal trilinear coupling A and tan (3 are uniquely fixed. Before investigating 
the impact of other values for the parameters m and A, let us first discuss the Higgs 
and sparticle spectrum as a function of M1/2 in the phenomenologically allowed region. 
The Higgs, neutralino and stau masses are shown in Fig. 2 where we also indicate the 
corresponding values of A . 

The essential features of the Higgs spectrum are as follows. For Mu 2 < 660 GeV, 
the lightest CP-even Higgs boson has a dominant singlet component, hence a very small 
coupling to the Z boson, which allows it to escape LEP constraints. The next-to-lightest 
CP-even scalar is SM-like, with a mass slightly above 114 GeV. The lightest CP-odd scalar 
is again singlet-like, with a mass above ~ 120 GeV. The heaviest CP-even and CP-odd 
scalars are practically degenerate with the charged Higgs boson, with masses above ~ 520 
GeV. For M1/2 ^ 660 GeV, the lightest CP-even scalar is SM-like with a mass increasing 
slightly with M\/2 up to ~ 120 GeV, while the next-to-lightest one is now singlet-like. 

-A (GeV) -A (GeV) 



100 150 200 250 300 350 400 100 150 200 250 300 350 400 




M 1/2 (GeV) M 1/2 (GeV) 

Figure 2: The Higgs (left) and neutralino plus stau (right) mass spectra in GeV as a 
function of M1/2 along the DM line; the values of A$ are indicated in the upper axis. 

The right-hand side of Fig. 2 displays the neutralino and stau spectrum, the lighter 
stau Ti being the next-to-LSP with a mass ~ 3 to 5 GeV above the one of the Xi singlino- 
like LSP. X2 an d X3 are ) respectively, bino and wino-like while the nearly degenerate x1,5 
states are Higgsino-like. The charginos Xi an d xf are nearly degenerate in mass with, 
respectively, X3 an d X45 (notice that the mass of the latter is ~ /i c ff)- The remaining 
sparticle spectrum is very "cMSSM"-like and can be obtained by running the program 
NMSPEC [9] with input parameters as in Fig. 3 below (and mo ~ 0) and also by using 
any cMSSM-based code, since the singlet sector practically decouples from the SUSY 
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spectrum. One approximately obtains m g 
generation) squark masses. 

Let us now discuss the impact of other values of the parameters m and A. As already 
stated above, the Higgs and sparticle spectra change very little with A provided that A 
remains small enough. Upper bounds on A result from LEP constraints on Higgs scalars 
with masses below the SM limit of ~ 114 GeV. For M\i 2 < 660 GeV, increasing A increases 
the mixing of the singlet-like CP-even scalar with doublet-like CP-even scalars and hence 
its couplings to the Z boson, which must not be too large. For My 2 ^ 660 GeV, a stronger 
mixing among the CP-even scalars can lower the mass of the lighter Higgs boson which is 
now SM-like, until it violates the LEP bound. 
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Figure 3: tan/5 and the upper bound on A as a function of M\/ 2 along the DM line. 



Fig. 3 shows the corresponding upper limits on the parameter A, which are particularly 
strong in the "crossover" -region near Mi/ 2 ~ 660 GeV, where relatively small values of A 
can generate a large mixing angle; in all cases, one has A < 0.02. For completeness, we 
also show the values of tan (3 along the DM line. 

We have also investigated the cNMSSM parameter space for non-zero ttiq and found re- 
gions which comply with all phenomenological constraints, including the DM relic density. 
However, in order to generate a non- vanishing vev of the singlet, the inequality mo < || A)| 
has to hold [2] and, thus, solutions only exist for small mo- Even for M\/ 2 ~ 1.5 TeV, 
i.e. for squark and gluino with masses of order 3 TeV and already beyond the LHC reach, 
only values for tuq below ps 140 GeV are allowed. Smaller M 1 / 2 corresponds to smaller A 
and, for instance, only < 20 GeV is viable for Mi/ 2 ~ 400 GeV. Since mo <C M\i 2 in 
all cases, we expect the phenomenology of the model to be similar to the one discussed 
here. More details will be given in a forthcoming publication [18]. 
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Finally, let us note that we have also verified possible instabilities of the potential 
along charge and colour breaking directions in field space, which could be dangerous for 
m « [2, 16]. Using analytic approximations to the potential along such dangerous 
directions, the inequality (A — ^Mi/2) < 9ml + §Mjy 2 was derived in Ref. [2], and 
is satisfied within the allowed region. Concerning directions in field space in which the 
potential is unbounded from below, it was shown in Ref. [17] that even if they exist, the 
decay rate of the standard vacuum is usually much larger than the age of the Universe. 

In conclusion, we have shown that the NMSSM with universal boundary conditions at 
the GUT scale is a very constrained scenario. For small values of the universal scalar mass 
m , which are theoretically interesting and excluded in the cMSSM, all present collider 
constraints on sparticle and Higgs masses are satisfied. Moreover, the requirement of a 
correct relic density for the dark matter candidate further constrains the parameter space 
to a one-dimensional [M1/2, Aq] line. Thus, only one single parameter, which can be taken 
as M1/2, is required to describe the salient features of the model. 

This model leads to an interesting phenomenology. For large Mi/ 2 , the lightest CP-even 
Higgs boson is SM-like with a mass smaller than ~ 120 GeV. For small M1/2, the lightest 
CP-even Higgs boson is a possibly very light singlet-like state which will, however, be very 
difficult to detect given the small value of A. In the SUSY sector, the singlino-like LSP 
will considerably modify the decay chains of sparticles [19]: one expects that all sparticles 
decay via the lightest tau slepton which then decays into the singlino-like LSP, leading 
to missing energy. For very small A, the lifetime of the tau slepton can become so large 
that its track can be visible [19]. In any case sparticle decays will differ in a spectacular 
way from the ones expected within MSSM-typical scenarios, which will allow to test the 
present scenario hopefully in the near future at the LHC. 
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